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CN ■ Abstract 
> 

We use the massless Thirring model to demonstrate a new approach 
' to non-perturbative fermion calculations based on the spherical field for- 

malism. The methods we present are free from the problems of fermion 
00 ' doubling and difficulties associated with integrating out massless fermions. 

' Using a non-perturbative regularization, we compute the two-point corre- 

lator and find agreement with the known analytic solution. 

n ' 

1 Introduction 

QJ I The massless Thirring model 0| is an exactly soluble system of a single self- 

i-C ' interacting massless fermion in two dimensions. There are a number of solu- 

tions in the literature based on properties of the Euler-Lagrange equations and 
fermion currents or bosonization techniques Given its simplicity and 

solubility, the model has become a popular testing ground for new ideas and 
■ methods in field theory. 

From a computational point of view, however, the massless Thirring model 
still presents a significant challenge. In the lattice field formalism, the difficul- 
ties are due to the appearance of fermion doubler states and singular inversion 
problems associated with integrating out massless fermions. In this work we 
use the model to illustrate new non-perturbative methods in the spherical field 
formalism 0-0. The techniques we present are quite general and can also be 
applied to other modal expansion methods such as periodic field theory |l0|| . 

As noted in |8|, we will not need to worry about fermion doubling. This 
is true for any modal field theory and follows from the fact that space is not 
discretized but retained as a continuous variable. Since our model is not super- 
renormalizable we will use a procedure called angle-smearing, a regularization 
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method designed for spherical field theory [gj . With angle-smearing regulariza- 
tion and a small set of local counterterms, we are able to remove all ultraviolet 
divergences in a manner such that the renormalized theory is finite and trans- 
lationally invariant. Comparison of our results with the known Thirring model 
solution will serve as a consistency check for our regularization and renormal- 
ization procedures. 

The organization of this paper is as follows. We begin with a short summary 
of the massless Thirring model, following the solution of Hagen Q|. Using angle- 
smearing regularization we obtain the spherical field Hamiltonian and construct 
a matrix representation of the Hamiltonian. We reduce the space of states using 
a two-parameter auxiliary cutoff procedure. In this reduced space we compute 
the time evolution of quantum states using a fourth-order Runge-Kutta-Fehlberg 
algorithm. We calculate the two point correlator for several values of the 
coupling and find agreement with the known analytic solution. 



2 The model 

We start with a list of our notational conventions. Our analysis will be in two- 
dimensional Euclidean space, and we use both cartesian and polar coordinates, 

t = {tcos9,tsm9) — {x,y). (1) 

The components of the spinors ip and ip are written as 

i, = [ iP'^ ] . (2) 

Our representation for the Dirac matrices is 



7 = i^, (3) 

and so 7 satisfies 

{f,7'}=-25'^ hj = h2. (4) 
The massless Thirring model is formally defined by the Lagrange density 
^ = #7-V^-|J-J, (5) 

where j is the fermion vector current. Johnson ||^ emphasized the importance 
of defining the regularized current precisely, and this was further clarified by the 
work of Sommerfield j^] and Hagen . We will use a regularization technique, 
introduced in called angle-smearing. We define the regularized current as 

/ = 5 i^s'ytps - Trljtpstps]) , (6) 
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where 



ijs{t,e) = ^ I de^{t,e + e). (7) 

We identify the radial variable t as our time parameter, and our definition of 
the current is local with respect to t. 

Hagen described the solution of the Thirring model in the Hamiltonian 
formalism with currents defined as products of the canonical operators at equal 
times. Though our equal-time surface is curved, the curvature of the integration 
segment in scales as while the ultraviolet divergences in this model are 
only logarithmic in M . In the M oo limit we therefore recover the standard 
results. As discussed in [Q, there exists a one parameter class of solutions to the 
Thirring model depending on the preferred definition of the regularized vector 
and axial vector currents. We will use the conventions used in and S, which 
in Hagen's notation corresponds with the parameter values ^ = i] = With 
this choice the Hamiltonian density takes the form 

n = Hfree + ■ jf + ■ J)^ (8) 

where 



27r ■ 



(9) 



The only counterterms needed in this model are wavefunction renormaliza- 
tion counterterms, a result of our careful definition for the regularized interac- 
tion in (^) . As in Q we calculate correlation functions using an unrenormalized 
Hamiltonian. The divergent wavefunction normalizations will appear simply 
as overall factors in the correlators. 



3 Spherical field Hamiltonian 



In this section we derive the form of the spherical field Hamiltonian. We first 
expand the fermion current in terms of components of the spinors, 



e-*' 



(10) 



-iS 



The anti-commutators of the regulated fields areg 



(11) 



(12) 



^Our definition of the Euclidean fermion fields and anti-commutation relations follows the 
conventions of [pj| . 
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{^i, V'J} = \ /ie*)de = A{t)e 



i0 



where 



(13) 
(14) 
(15) 



A(i) = ^sin(^). 
From, the anti-commutation relations, the t component of the current is 

and so 

{i ■ jf = A{t) [e-^"^l^i + e'ViVl] - ZVJIV'iVii^J - A^{t). (16) 
Similarly we find 

{e ■ j f = A{t) [e-^Vlvi + e^VivI] - z^lvi^iV'l- (17) 

The Hamiltonian can now be written as 

H = Hfree + J d0 1 [A(i) [e'^^^l^i + e''^ii>l] - 2i,li,ii,i^l] } . (18) 

We have omitted the constant term proportional to A'^{t). 
Let us define the partial wave modes 



(19) 
(20) 



It is straightforward to show that for n ^ 0, 



\Mt) \Mt ) 

We can extend this result to the case n = using the convenient shorthand 



Mi)- 



(21) 



\ Mt ) 

At this point it is convenient to rescale "0, 
In terms of the partial waves, 



n 

-iE 



(22) 
(23) 
(24) 



E 

-ni+7i2— n3+n4=0 



("Mt)(T^T 

jV^l'„,V'i+iV'-„3-iV'. 



V Mt A Mt J\ Mt }\Mt} 
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where 



b = T^- (25) 

Since b is the parameter appearing in the Hamiltonian, it is somewhat more 
convenient to express c in terms of 6, 



^ ~ b 

Let us define the ladder operators^ 
These operators satisfy the anti-commutation relations 



(26) 

(27) 
(28) 



(29) 



with all other anti-commutators equal to zero. We can now recast the Hamil- 
tonian as 



H 



bTTtAjt) sin(-jg^) sm(^) 
(mT ) (~mt) 



(30) 



E 

-ni+n2+n3— n4— 



2 "n3 "n4 (- "l-i y "2 U "3-i V 

V A'/t A Mt / V Mt J\ Mt J 



We will implement a high spin cutoff by removing terms in the interaction 
containing operators a^, a]j, a^^ , or aj^''' for |n| > Jmax- This has the effect of 
removing high spin modes, which correspond with large tangential momentum 
states. We should emphasize, however, that Jmax is an auxiliary cutoff. It does 
not play a role in the regularization scheme since the interactions have already 
been rendered finite using angle-smearing. The contribution of these high spin 
modes is negligible so long as 



> M, 



(31) 



where t is the characteristic radius of the process being measured. Returning 
back to and (|l|) and removing the contribution of these partial waves, we 
find that A{t) is replaced by 



A, 



At) 



Jmax 
^ V 

27rt / J 
"i\— — J 



( ITf ) ( Mt ) 



(32) 



This notation is slightly different from that used in 



The translation is as follows; 



J. it _ i- i+. T Tt _ T- T 
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Let \0) free ground state of the free massless fermion Hainiltonian.| 

For n > 0, we find 



4 \0)free = 4 \0) free = (33) 

and for n < 0, 

4^|0)/„e=4^|0)/.e.=0. (34) 

It is convenient to define the normal-ordered products 

it i _ / "-n^-n for n > Tt T - / "'n^'i for ?^ > , . 

- \ -aia^J torn < " \ -alaTt for n < 0. ^"^^^ 

The ordering for other operators is immaterial since the anti-commutators are 
zero. We can now rewrite H in terms of normal-ordered products, 

- (f + 0{J,~L)) (ai^ai + aTtaT ) (36) 



-rii+n2+"3-"4=0 '" 



V Mt }\Mt}\ Mt )\Mt} 



There is an 0{J^^^ term due to a small asymmetry in our cutoff procedure 
with respect to the two boundaries — Jmax and Jmax-^ We will neglect this term 
in the limit Jmax ^ c>o- 

4 Two-point correlator 

We wish to study the properties of the two-point correlator. The massless 
Thirring model is invariant under the discrete transformation 

^i(t),^T(i)^_^i(i)^_^T(i)^ (37) 
as well as the transformation 

^l;Ht,0),i>\t,d) ^ i^\t,-e),i^'^{t,-d). (38) 
From these we deduce 

(0| T |0) = (0| T [V^HilV'HO)] |0) = (39) 



•^The ground state of the free massless Hamiltonian is actually degenerate due to s-wave 
excitations, but this is remedied by taking the m > limit of the massive fermion theory. 

^If desired we could eliminate this term and the asymmetry by a slight change in the 
angle-smearing procedure for {[). 
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and 



{0\T[^'^{t,e)ip^iO)] \0) ^ {0\T [il;-^{t, -9)i;^ (0)] |0) . (40) 

It therefore suffices to consider just the correlator on the left side of (|40|). 
In the limit M ^ oo the form of the correlator is given by 

{0\T[^^{t,9)^^{0)] |0) = ^{k{c)M)f^t^, (41) 

where k{c) is a dimensionless parameter. Standard analytic methods do not 
yield a simple closed form expression for k{c). We will therefore extract fc(c) 
from the computed value of the correlator at a specific renormalization point 

We define 

|0). (42) 
Since 



f{t) = {0\T 



{O\T[^Ht,0)ij^{Q)]\O)^^{O\T i^l{t)^U0) |0), 



(43) 



we conclude that 



fit) = {k(c)M)^t 1- 



(44) 



We now compute f{t) using the spherical field Hamiltonian. We first need a 
matrix representation for the Grassmann ladder operators. We will use tensor 
products of the 2x2 identity matrix and Pauli matrices: 



4 = cr^ 0-2 «) {\(T^ + |cr.y) 1, (45) 

4 = ® ^z® {\(Jx + f CTy) 10 1. 

The representations for a]^ and a]^ are defined by the conjugate transposes of 
these matrices. We can now calculate the correlator J[t) using the relation [p| 



f(t) = lim lim ^ ' r X TTTr^ ^_J_JiL_ UQ) 



A straightforward calculation of (fig), however, is rather inefficient. There are 
several techniques which we will first use to simplify the calculation. 

^In some regularization schemes fc(c) can be calculated analytically ^[^, and it may be 
worthwhile to use these techniques in future work. In this first analysis, however, we prefer 
to present a more straightforward and typical example of the angle-smearing regularization 
method. 
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The time evolution of the system at large t is dominated by the contribution 
of the ground state or, more precisely, the adiabatic flow of the t-dependent 
ground state. As discussed in |^ a similar phenomenon occurs at small t, due 
to the divergence of energy levels near t = 0. It is therefore not necessary to 
compute the full matrix trace in the numerator and denominator of (^6|). It is 
instead sufficient to compute the corresponding ratio for a single matrix element. 
After making this reduction, we can then go a step further and eliminate states 
which do not contribute to the matrix element. 

The high spin parameter Jmax was used to remove high spin modes with 
\n\ > Jmax- This, however, is not a uniform cutoff in the space of states and 
most of the remaining states are still high kinetic energy states. Although none 
of the individual modes are energetic, many of the modes can be simultaneously 
excited. Let us define and A^^ to be bit switches, 1 or 0, depending on 
whether or not the corresponding mode is excited. Let us also define a cutoff 
parameter Kmax- We will remove all high kinetic energy states such that 

E„ [H (A^i + ND] > i^max. (47) 

For consistency A'max should be about the same size as Jmax- 

5 Results 

The CPU time and memory requirement for calculating (|4^) scales linearly with 
the number of transitions in H (i.e., non-zero elements in our matrix represen- 
tation). In Table 1 we have shown the number of states and transitions for 
different values of Jmax- 



Tablc 1 



Jmax 


2 


4 


6 


8 


10 


12 


states 


6 


40 


210 


920 


3600 


13000 


transitions 


38 


500 


4200 


26000 


1.4E5 


3.9E5 



We have calculated f{t) for Jmax < 12 and several values of the coupling b. The 
total run time was about 100 hours on a 350 MHz PC with 256 MB RAM. 

The matrix time evolution equations in (^) were computed using a fourth- 
order Runge-Kutta-Fehlberg algorithm. We have set 

^max — 2 Jmax + 2. (48) 

We will use the notation fj^^^ (t) to identify the corresponding result for a given 
value of Jmax- In Figure 1 we have plotted f.j^.^^{t) for 6 = 1 and Jmax — 4, 8, 12. 
We have scaled t and /(<) in dimensional units chosen such that M = 3. For 
finite Jmax we expect deviations from the Jmax oo limit to be of size 0{ ^j.i * ). 
The curves shown in Figure 1 appear consistent with this rate of convergence. 
We can extrapolate to the limit Jmax oo using the asymptotic form 

/J„„. it) = /oo it) + Jma'x/''^ (*) + ■••• (49) 



8 



For 6 = 0, 0.5, 1.0, 3.0 and M — 3 we have calculated f{t) using this extrapola- 
tion technique for Jmax = 10 and 12]^ The results are shown in Figure 2. For 
comparison we have plotted the analytic solution 

/^(i) = (fc(c)M)T^t^. (50) 

The relation between b and c can be found in (|2^) and (p6|). The parameter 
fc(c) is fitted to the value of the correlator at the renormalization point t = 0.6.|] 
The agreement appears quite good. Some deviations from the analytic solution 
are due to 0{j^^) residual terms, which were left out of the derivation of (pO|). 
These effects are significant in the small t region, t < M~^, especially for larger 
values of b. The values we find for fc(c) are shown in Table 2.^ 



Tabic 2 



b 


0.5 


1.0 


3.0 


He) 


1.77 


1.77 


1.68 



We can compare our results at small b with a simple perturbative calculation. 
Evaluating the corresponding regulated two-loop diagram we obtain, for small 
b, k{c) « 1.75. This appears consistent with the results in Table 2. 



6 Summary 

We derived the angle-smeared spherical Hamiltonian for the massless Thirring 
model and constructed an explicit matrix representation. We discarded negligi- 
ble high energy states using auxiliary cutoff parameters Jmax and if max- In this 
reduced space we computed the time evolution of quantum states and calculated 
the two-point correlator for several values of the coupling. The results of our 
computation are in close agreement with the known analytic solution. In addi- 
tion to demonstrating new computational methods, our analysis also serves as a 
consistency check of the regularization and renormalization methods introduced 




We believe that this work represents a significant new direction in the non- 
perturbative computation of fermion dynamics. Future work will study the 
application of these methods to systems of interacting bosons and fermions. 
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^Both our results and the analytic solution are even in b, and so we consider only positive 
values. 

^The relative error is expected to be small in the vicinity of this point. 

*In some regularization schemes fc(c) can be shown to be independent of the coupling. Our 
regularization method seems to be rather close to this, with only a slow variation with respect 
to the coupling strength. 
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Figures 

Figure 1. Plot of /j_(i) for 6 = 1 and J^ax = 4, 8, 12. 
Figure 2. Plot of /^(t) and f(t) for 6 = 0, .5, 1, 3 and M = 3. 
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Figure 1 
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